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ABSTRACT 

Quantitative models of the external magneto spheric field have been 
derived by making least^squares fits to magnetic field measurements from 
four IMP satellites. The data set consists of 12,616 vector field averages 
over half-R^ intervals between 4 and 17 taken from 451 satellite 

orbits between 1966 and 1972. The data were fit to a power series expan- 
sion in the solar magnetic coordinates and the solar wind-dipole tilt 
angle, and thus the models contain the effects of seasonal north-south 
asymmetries. The expansions are divergence-free, but unlike the usual 
scalar potential expansions, the models contain a non-zero curl representing 
currents distributed within the magnetosphere. Characteristics of four 
models are presented, representing different degrees of magnetic dis- 
turbance as determined by the range of Kp values. The latitude at the 
earth separating open polar cap field lines from field lines closing 
on the day side is about 5^ lower than that determined by previous 
theoretically-derived models. At times of high Kp, additional high- 
latitude field lines are drawn back into the tail. Near solstice, the 
separation latitude can be as low as 75^ in the winter hemisphere. 

The average northward component of the external field is much smaller 
than that predicted by theoretical models, indicating the important 

effects of distributed currents in the magnetosphere. Current densities 

-9 2 

implied by the models are of the order of 10 amp/m across the magnetotail. 
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INTRODUCTION 

Quantitative models of the geomagnetic field generally fall into 
two types: those that describe the field due to sources inside the earth 

(internal models) and those which attempt to describe the effect of currents 
located above the earth’s surface (external models). Internal models are 
derived from direct measurements of one or more components of the field. 

Since the time of Gauss, spherical, harmonic expansions of a scalar potential 
have been used to fit these measurements. Similar techniques have been used 
to model the overhead current systems responsible for the diurnal and lunar 
variations in the surface field. However, up until now quantitative models 
of the external magnetospheric field have not been based on field measure- 
ments or are based only indirectly on them. Some of these models ( Mead , 

1964; Midgley , 1964; Olson , private communication; Choe and Beard , 1973) 
have been based on theoretical solutions to the traditional Chapman-Ferraro 
problem of an unmagnetized plasma incident upon a magnetic dipole (e.g* , 

Mead and Beard , 1964; Olson , 1969; Choe et al . , 1973)- In these models the 
resulting external field is expressed in terms of a scalar potential, 
expanded in spherical harmonics and modified to reflect the fact that the 
current sources are external to the region of the expansion. Others ( Hones , 
1963; Taylor and Hones , 1965; Williams and Mead , 1965; Antonova and 
Shabanskiy , 1968; Sugiura and Poros , 1973; Olson and Pfitzer , 1973) have 
been based partly on theory and partly on known characteristics of the 
external field. These models have -used some combination of spherical 
harmonics, image dipoles and/or current sheets to describe the magnetospheric 
field. The characteristics of many of these models have been discussed by 


Roederer (1969) . 


preceding page blank not filmed 
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In this study we have made two major departures from previous work. 
First of all, we follow procedures heretofore used only in modelling the 
internal field, namely, making least- squares fits to a body of magnetometer 
data using a model with adjustable coefficients. The data set consists of 
magnetometer measurements from four spacecraft; Explorer 33 (A-Imp-D) , 
Explorer 34 (Imp 4), Explorer 41 (Imp 5), and Explorer 43 (Imp 6 ), These 
spacecraft were in highly eccentric orbits, with apogee ranging from 29 
to 80 Rg. The data were obtained over an interval of almost six years on 
451 separate spacecraft orbits. 

Secondly, we found it necessary to abandon the use of a scalar 
potential with its implicit assumption that there are no current sources in 
the region of the expansion. Instead, we used power series expansions of 
the field components themselves in the solar magnetic coordinates and the 
tilt angle. Although VxS 7 ^ 0 in our expansions, indicating the presence of 
distributed currents, the fitting program constrained the coefficients so 
as to make = 0 . 

The coefficients and characteristics of four models are presented in 
this study. These models contain terms up to quadratic in the solar magnetic 
coordinates and linear in tilt angle. Each was derived by least-squares fits 
to subsets of the data sorted according to the Kp value- The models contain 
only terms which retain the obvious east-west and summer -win ter symmetries. 

By taking the curl of the resulting expressions for the magnetic 
field, a general picture of the currents can be determined. Although the 
models do not contain sufficiently high-order terms to determine these 

_9 2 

currents with much precision, current densities of the order of 10 amp/m 
in the outer magnetosphere are implied by the models. 
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since the models are derived from measurements of the external 
field, the accuracy of the models is limited by the spatial coverage of 
the data. No data was used beyond 17 R^, and therefore the models should 
not be used to calculate magnetic fields beyond this distance. Very 
little data was available inside 5 Rg, and therefore the models do not 
accurately represent the field depressions observed by Sugiura (1973) and 
others at 2-5 R^. In addition, data were often not available at some 
longitudes during certain seasons. These gaps in the data coverage, 
together with the constant variability in magnetospheric configuration, 
limited us to very simple models - 

In this paper a number of model characteristics are compared with 
experimental results previously presented by Fairfield (1968, 1971) , 
Sugiura et al . (1971), Sugiura (1972a), and Skillman (1973). In a 
companion paper ( Fairfield and Mead , 1973) the models are used to map 
field lines between low altitudes and the distant magnetosphere so as to 
compare, for example, the observed positions of the high*-latitude 
boundary of the trapping region with the polar cap field lines as 
determined by the various models. 
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TREATMENT OF DATA 


The magnetic field data used as the starting point in this study 
Were averages of the three vector components over telemetry seguences 
of duration 81 sec (Explorer 33), 20.5 sec (Explorers 34 and 41), and 
15.4 sec (Explorer 43). These experiments have been described previously 
( Behannon , 1968; Fairfield , 1969; Fairfield and Ness , 1972; Fairfie_ld , 
1973a) ; all were performed on magnetically clean spacecraft, yielding 
vector field measurements with an absolute accuracy of better than 
1 gamma. The least sensitive ranges of the four instruments, which 
determined the earthward termination of the data, were ±64y (Explorer 33) , 
±128t (Exlorer 34), ±200y (Explorer 41), and ±432y (Explorer 43). Thus 
no measurements were available inside 4 and relatively few inside 5 

Plots of both the measured sequence averages and the theoretical 
predictions from a model of the internal field were used to identify 
magnetopause crossings and to detect and eliminate measurements which were 
contaminated by telemetry noise. (For examples of these plots see 
Fairfield and Ness , 1972; Fairfield , 1973b.) Data taken outside the 
magnetopause were eliminated for the present analysis. 

Since both time and spatial variations were small between sequence 
averaging intervals, further averaging was clearly desirable in order to 
minimize redundant data and reduce the data set to a manageable size. In 
order to minimize the effect of variations in the internal field over the 
averaging interval, the quantities AB, AI, and AD were calculated for 
each telemetry sequence, where Ab is the difference in magnitude between 
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the measured and internal reference field, and AI and AD are the inclina- 
tion and declination differences ( Mead and Cahill , 1967) . The internal 
reference field was, in most cases, the IGRF model ( lAGA ^ 1969). These 
differences were then averaged over half-earth-radii intervals of radial 
distance. Typical averages were over 10 to 15 minutes of time, thus 
suppressing the most rapid time variations. (The orbit of Explorer 33 
was more circular than the others, and there we also required that the 
averaging interval not exceed 5^ in longitude.) Associated with each of 
these data points was the average spacecraft position in solar magnetic 
coordinates, the tilt angle T, and the time of the measurement. Only 
data within 17 was included in this analysis, as we felt that our 
simple model expansions would have difficulty in simultaneously fitting 
the near-earth data and the distant tail field data. 

The extent of the data set is described in Table 1 . The data were 

obtained over an interval of almost six years, covering 4340 spacecraft 

hours within the magnetosphere on 451 separate spacecraft orbits. Over 
7 . 

4 X 10 individual vector measurements were combined into 12,616 half- 
earth- radii averages. 

The data were then separated into two roughly equal parts, depending 
on the Kp value at the time each measurement was made. The quiet data 
included measurements for which Kp < 2; for the disturbed data, Kp > 2. 
Subsequently, two more limited subsets were created, a super-quiet (Kp = 0, 0^ 
data set and a super-disturbed (Kp > 3) set- 
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To prevent regions with a high concentration of data from having 
a disproportionately large effect on the models/ we decided to carry out 
further averages over our four-dimensional space* Intervals of solar 
magnetic latitude and longitude were chosen so as to divide a sphere 
around the earth into 96 equal solid angles (see Fairfield and Ness / 1967, 
for definition of the angles) . This division, together with 13 radial 
intervals of 1 (4 < R < 17 R^) and 7 tilt-angle inteirvals of 10° 

(-35° < T < 35°) defined 96 x 13 x 7 -= 8736 four-dimensional "buckets" - 
Each data point was then sorted into the appropriate bucket and the 
spatial parameters and field values within a given bucket were averaged. 
Four final data sets resulted: the quiet-time data, with 2368 points; 
the disturbed data, with 2206 points; the super-quiet data, with 882 points 
and the super- disturbed data, with 1334 points* 

This data, however, was by no means uniformly distributed over our 
four-dimensional space- Each spacecraft traverses a limited range of 
latitudes and, furthermore, makes measurements within a given longitude 
sector only during certain times of the year, when the range of tilt angles 
is correspondingly limited. The resulting non-uniform distribution of data 
can only be minimized by having a large number of spacecraft with a variety 
of orbital configurations. 

If the quiet data had been evenly distributed, about 85% of the 
buckets inside the magnetopause would have contained a single point. In 
fact, however, only 31% of these buckets were filled, with a maximum of 
12 points and an average density of 2.7 points per bucket. The empty 
buckets are at latitudes not sampled by the various spacecraft and at 


certain longitudes not san^led for certain tilt angles. 
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The distribution of data is illustrated in Figure 1 ^ which indicates 
the magnetic latitude and local time of each point in the final data set 
(quiet and disturbed combined) for distances of 6-8 (top) and 13-15 
(bottom) - The near-earth data indicates a scarcity of points near the 
magnetic equatorial plane and better latitudinal coverage in the northern 
hemisphere compared to the southern hemisphere • There are practically no 
measurements south of -50^ latitude- At the greater distances, data is 
absent near the noon meridian because such measurements would be outside 
the magnetopause- Fewer high-latitude data are present at greater 
distances, because the spacecraft are closer to apogee, which for all 
spacecraft is in the vicinity of the solar ecliptic equatorial plane- 
The various symbols indicate tilt angles within three different 
ranges- It is clear that coverage in certain latitude-local time regions 
is restricted to certain seasons- Particularly apparent in the near- 
earth distribution is a lack of southern hemisphere data in the nightside 
hemisphere for positive tilt angles and in the dayside hemisphere for 
negative tilt angles- These gaps in the four-dimensional data distribution 
mean that a least-squares model will be poorly constrained in regions 
where there is no data- This situation imposes a severe limitation on 
these models, a point which will be discussed further in a later section- 
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RESULTS 


In the initial phase of the study we attempted to fit the measured 
difference field with spherical harmonic expansions of the scalar potential 
similar to those used by Mead (1964) in his theoretically-derived model of 
the external field. The values of the coefficients resulting from fits to 
different components r however, were very different from each other. We 
concluded that the assumption that the region was source-free was simply 
not valid; distributed currents within the region of measurements, which 
cannot be represented by a scalar potential, made a substantial contribu- 
tion to the difference field. 

Since Mead (1964) had noted that the expansion of the external 
field was particularly simple when expressed in cartesian coordinates , we 
transformed the difference field into cartesian coordinates (AB^, AB^, and 
AB ) and modified our least-squares program so as to fit each component 
to a power series expansion in cartesian solar magnetic coordinates, 
including terms linear in the tilt angle. (In this coordinate system the 
z-axis is aligned with the magnetic dipole, the x-z plane contains the 
earth-sun line, positive values of x are on the dayside and the y-axis is 
in the dusk meridian. The tilt angle, T, is the complement of the angle 
between the z-axis and the earth-sun line, or, equivalently, the geomag- 
netic latitude of the subsolar point; it is positive during northern 
hemisphere summer.) 

Neither 7‘B nor 7 xb necessarily vanished with this expansion. 
However, the coefficients obtained for a quadratic fit (20 coefficients 
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for each component, totaling 60) were found to come very close to satis- 
fying thev*^ = 0 condition, although Vx§ was clearly non-zero. This low 
-y 

value ofV*B gave us great confidence in the self-consistency of the under- 
lying data (although since the data was not taken at one instant of time 
there was no reason to expect V-S to be identically zero) . 

Thus encouraged r we developed a more general least-squares program/ 
using a method of Lagrangian multipliers/ to fit all three components of 
the field simultaneously, subject to the condition that the coefficients 
satisfy the 7«B 0 requirement , (For the quadratic fit, eight equations 

of constraint were imposed; thus there were only 52 independent coeffi- 
cients. ) The imposition of the V*B = 0 condition changed the values of 
the largest coefficients by only a few percent and the residuals increased 
by less than 1% (see below for definition of the root-mean-square vector 
residual field) . 

Two 60-coefficient models resulting from this program have been 
made available to a number of workers and have been used as the basis for 
a few subsequent papers. In these models all possible terms through 
quadratic in the solar magnetic coordinates and linear in the tilt angle 
were included/ with no symmetry imposed. The models/ however/ showed a 
high degree of both dawn-dusk and north-south symmetry. The dawn-dusk 
symmetry was greatest about a meridian plane rotated about 4^ from the 
solar direction/ an angle about equal to the aberration of the solar wind 
direction caused by the orbital motion of the earth. Also, the field-line 
topology in one hemisphere for a given positive tilt angle was very similar 
to the topology in the opposite hemisphere for the same negative tilt. It 



12 


was not possible, however, to determine whether the small dawn-dusk and 
north-south asymmetries represented real asymmetries in the magnetosphere, 
or whether they were caused by time variations and non-uniformities in 
the distribution of data. Since we had no assurance of the reality of 
the model asymmetries, we developed a modified least-squares program to 
fit an expansion containing only those terms which retained both types of 


symmetries. The dawn-dusk symmetry condition requires that 

ABjj(-Y) = AB^(Y) (1) 

ABy(-Y) = -ABy(Y) (2) 

AB^C-Y) = Ab^(Y) (3) 

for all X, Y, Z, and T. North-south symmetry requires that 

ABjj(-T,-Z) = -ABj^(T,Z) (4) 

ABy(-T,-Z) = -ABy(T,Z) (5) 

AB (-T,-Z) = AB,(T,Z) (6) 


These conditions (which are automatically satisfied for a dipole field) 
eliminate almost three-fourths of the coefficients, and the quadratic 
expansion becomes 


AB = 

2 2 2 

a^Z + S2XZ + T(a3 + a4X + a^X + a^Y + a^Z ) 

(7) 

AB == 

b^YZ + T(b2Y + b3XY) 

(8) 

= 

c^ + C 2 X + + c^Y^ + C5Z^ + T(CgZ + CyXZ) 

(9) 


(note comments on units below), subject to the V*B — 0 conditions: 

a_ + b, + 2Cc = 0 (10) 

a^ + b2 + Cg = 0 (11) 

2a^ + b^ + Cy = 0 (12) 

This is a 17-coefficient expansion, with three restrictions, making 
14 independent coefficients. 
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Since the preliminary w6rk had indicated that the optimum plane of 
symmetry would be rotated a few degrees from the noon-midnight meridian# 
the X and Y coordinates and field components were rotated through an 
arbitrary angle before determining the least-squares coefficients to the 
symmetric expansion • Residuals were then calculated for various rotation 
angles, A minimum in the residuals was reached at a rotation angle of 
from 3 to 5 degrees# depending upon the data set# beyond which the 
residuals began to increase. For the more disturbed data sets# the 
rotation angle for minimum residuals was less than for the quiet data 
sets, suggesting that the aberration angle is less during disturbed times, 
when the average solar wind velocity is probably higher. Our results 
do not support the suggestion made by Cummings et al , (1971) # based on 
analysis of ATS-1 magnetometer data, that the plane of maximum symmetry 
is rotated by 10-15^- We found the residuals for such large rotation 
angles to be significantly larger than for angles of 3-5^. 

The residuals for the 17-coefficient symmetric model with a 
rotation of 4° were only about 4% larger than the residuals for the 
60-coefficient model with no symmetry imposed. The fact that elimination 
of over seventy percent of the coefficients only increased the residuals 
by this small amount gave us added confidence in the basic symmetry of 
the underlying data. 

In addition to a quiet and a disturbed model# coefficients were also 
determined for a super-quiet and a super-disturbed model. The coefficients 
(in gammas) for all four of these models are listed in Table 2# valid for 
X, Y# and Z given in units of tens of earth radii (1 = 6378 km) and T in 

units of tens of degrees. These units were chosen so as to make the 
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contribution from each term, in gammas r roughly proportional to the 
tnagnitude of the corresponding coefficient for moderate tilt angles (-10^) 
at geocentric distances of about 10^15 R^. 

Perhaps the most obvious feature of the coefficients is the tendency 
for them to increase in magnitude with increasing Kp, indicating generally 
larger magneto spheric distortions during times of enhanced magnetic activity. 
It is particularly interesting that the c^ coefficient has a negative sign. 
This term is equivalent to the g^ term in the Mead (1964) model, which gave 
a northward field of 25y for a sub-solar boundary distance =10 Rg. 

The Cj_ coefficient in these models, however, gives a southward field of 
10-197, depending on the Kp range. Since currents at the magnetopause 
give rise to a northward-directed field, this result implies that the 
distributed currents within the magnetosphere give rise to a southward 
field near the earth which more than offsets the effect of the magnetopause 
currents. (The c^ and terms increase the AB^ component in the outer 
magnetosphere; thus the average measured value of AB^ was actually -1.6 y 
for the quiet data and -5.5y for the disturbed data.) The total field 
magnitude in the equatorial plane near the earth calculated from these 
models is much less than the average magnitude predicted by a model con- 
taining boundary currents only- This is consistent with Sugiura * s (1972a) 
results indicating negative values of AB near the magnetic equatorial plane 
close to the earth. It is clear that any realistic model of the magneto- 
spheric field must have the ability to describe distributed currents. 

The coefficient with the largest magnitude is a^, which gives a 
field directed generally towards the sun in the northern hemisphere and 
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away from the sun in the southern hemisphere r thus producing oppositely- 

directed fields on either side of the plasma sheet in the geomagnetic 

tail. Another large term is c^/ giving a northward field on the dayside 

and a southward field on the nightside. A restriction that the models 

be curl^free would have required = c^r and these coefficients would 

then be the equivalent of /J g^ - 21 y/ 10 R in the Mead (1964) model. 

Thus, the values of and in the MF73Q model agree quite well with 

.1 

Mead’s earlier value of g^/ and the fact that a^ > c^ is an indication of 
a magneto spheric current directed from dawn to dusk (see discussion of 
curl B in a later section) . 

The root-mean-square values of the measured vector difference 
field, vector residual field for each of these 

models is shown in Table 3 . These quantities are defined by: 

' Sint)^/Nl ’’ (13) 

1=1 

l“lres = ij^<4,eas " Bl„t “ 

where is the IGRF 1965.0 model of the internal field, ^^jj^odel 

given by Equations 7-9, and the quantities within the parentheses are 
vector differences, not differences between scalar magnitudes, as might 
be given if only scalar measurements were available. The summations 
are over the N bucket averages for each subset of the data. 

As might be expected, increases with increasing magnetic 

activity, the measured difference field being about twice as great for 
periods when Kp > 3 as when Kp = 0 or 0^. The residuals calculated with 
the appropriate models increase more or less proportionally. A figure of 


merit. 


defined by 
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F = AB / AB (15) 

' 'res ' 'rms 

which is a rough measure of the goodness of fit^ is relatively constant 
at about 0-5. Thusr only about half of the measured vector difference 
field can be predicted by the appropriate model- We believe that time 
variations in the magneto spheric field are responsible for most of the 
residuals. Since dynamic changes are continually taking place within 
the magnetosphere r measurements made at different times at the same 
spatial position with the same tilt angle and Kp value would most likely 
yield values of the vector field differing among themselves by amounts 
roughly comparable to the residuals given in Table 3, Thus, it is our 
opinion that even if reliable measurements were available throughout the 
entire magnetosphere during all seasons of the year, and even if a model 
were used containing many more terms than these models, the figure-of- 
merit F could not be reduced much below the values given in Table 3- A 
model of the internal field such as that of Cain and Sweeney (1970) fits 
the scalar field in the altitude range 400-1500 km (where |b| = 20^000- 
40,000y) with an rms scalar residual of 9y- It would appear that time 
variations throughout the magnetosphere, even during quiet times, are of 
the order of lOy^ and this represents about the limiting accuracy of any 
model which attempts to represent the average field configuration within 
the magnetosphere - 

It is instructive to review the effect which changes in the form 
of the model have upon the root-mean- square residuals- We take the 17- 
coefficient MF73Q model as a base with which to compare the residuals 
obtained with other possible models- This model includes all terms 
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through quadratic in the spatial variables and linear in the tilt angle 
which satisfy the symmetry requirements of Equations 1-6. Restricting 
this model to only linear terms (7 coefficients) increased the residuals 
by 24%. A 60-coefficient quadratic model with both symmetric and 
asymmetric terms gave residuals 4% lower than the symmetric model ♦ A 
30-coefficient quadratic model without any dependence on tilt angle 
gave residuals 30% higher than the 17-coefficient model. A 33-coefficient 
model which permitted all symmetric terms through cubic in the spatial 
variables and linear in the tilt angle reduced the residuals by 6%. How- 
ever, this cubic model exhibited a number of very unphysical characteristics 
in regions of 4-dimensional (X, Y, Z, T) space for which no data was 
available. 

Thus the 17-coefficient symmetric quadratic model, with coefficients 
restricted so as to satisfy the V*B = 0 condition, has the best overall 
characteristics of any of the models we have derived to date. It is quite 
clear that any accurate model of the magneto spheric field must include a 
dependence upon tilt angle and must reflect in some fashion the presence 
of distributed currents. It seems that a quadratic model is the highest- 
order realistic model that can be derived from the present data set. 

Before higher-order models are feasible, additional data coverage is 
required: more measurements at very high latitudes, more complete 

distributions of data in longitude at all tilt angles, and vector 
measurements over a wide latitude range at distances less than 4 R^. 
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MODEL CHARACTERISTICS 

Representative field lines (every 2° in latitude) for each of the 
four models in the noon-midnight meridian (more specifically, in the 
plane of symmetry rotated 4*^ from the noon-midnight meridian) for tilt 
angle T = 0° are shown in Figure 2 . Contours of constant total field 
magnitude are shown as dashed lines. There is a marked tendency during 
periods of high Kp for the equatorial crossing distance of the last 
closed field line on the noon meridian to be closer to the earth and for 
additional high-latitude lines to be swept back into the tail. In the 
super-quiet (SQ) model the 76*^ line crosses the equator at 11.3 R^ and 
the 80° line at 12.1 R^ on the dayside. In the super-disturbed (SD) 
model, the 76° line crosses at 10.4 Rg and the 78° and 80° lines are 
carried back into the tail. Thus, it appears that one effect of magnetic 
disturbance is to transport magnetic flux from the dayside back into the 
tail. The dayside latitude at the earth's surface separating the closed 
field lines from the polar cap lines is generally about 5° lower than that 
found in previous theoretically derived quantitative models, but is in 
good agreement with the experimental results of Fairfield (1968) . 

In each model there are two high- latitude neutral points in the 
noon meridian which separate those field lines which remain on the dayside 
from those that are swept back into the tail. These neutral points move 
closer to the earth and towards lower latitudes during periods of greater 
magnetic disturbance. In the SQ model they are located at = 62°, 

R = 13.5 Rg, and in the SD model at = 50°, R = 11.3 Rg. This compares 


o 
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with a neutral point latitude of about 72° in the theoretical model of 
Mead (1964) and 74° in the Olson (1969) and Choe et al - (1973) models, 

independent of solar wind intensity. 

An unexpected characteristic of each of the models is a minimum 
in the field magnitude (<10Y) at around 12“15 along the midnight 
equator, beyond which the field strength begins to increase. We believe 
that this is a non-physical characteristic, caused by the sparsity 
of magnetometer data near the equatorial plane, together with the fact 
that the expansions are only quadratic in the spatial variables and 
contain no explicit terms which would represent the effects of a current 
sheet in the tail. 

Figure 3 shows the noon-midnight topology and field magnitude 
contours for a tilt angle of 30°, corresponding to northern hemisphere 
summer. Here again there is a strong tendency during periods of enhanced 
magnetic activity for the last closed field line to be closer to the 
earth and for more field lines to be swept back into the tail. As is 
the case with almost all other tilted magnetosphere models, some of the 
field lines appear to penetrate the day side magnetopause, as though they 
were connected with the interplanetary field. 

One very clear effect is the change in latitude of the conjugate 
point which takes place near solstice. The dayside field line emerging 
at +76^ for T = 30°, for example, re-enters the earth at a latitude of 
-73-3^ in the SQ model and -71,4^ in the SD model. Thus, the dayside 
conjugate point in the winter hemisphere moves equatorward, sometimes 
as much as 500-600 km. In the MF73D model the 78^ line in the summer 



20 


o o , . 

hemisphere closes on the dayside, whereas both the 76 and 78 lines in 
the winter hemisphere are swept into the tail, and the 74^ line appears 
to penetrate the magnetopause* 

The effect is not as great on the nightside* For T - 30° the 
southern (winter) points conjugate to the +66° and +68° lines move about 
one or two degrees poleward in the various models* Field lines originating 
at 70° latitude and above trace beyond 20 in the tail, where the models 
are no longer applicaible* Conjugate point variations at the synchronous 
orbit are discussed by Fairfield and Mead (1973) * 

A broad minimum in the field magnitude (!b| < lOy) is located about 
midway between the solar magnetic and solar magnetospheric equators in the 
midnight meridian- There is again a (non-physical) tendency for the field 
magnitude to increase in value beyond 13-15 in the midnight equatorial 
region- The neutral points are retained in the noon meridian, but the 
winter hemisphere neutral points move towards lower latitudes. In the SD 
model the winter hemisphere neutral point (for T = 30°) moves from 50^ 
down to 41^ magnetic latitude, and the summer hemisphere neutral point 
moves up to 56*^ and inwards to 10-6 R^. These latitudes are substantially 
lower than those obtained by Olson (1969) or Choe et al . (1973) . The 
winter hemisphere neutral point in the model of Choe et al . moved from 
74° down to 64° for T = 30^; Olson found essentially no change in the 
magnetic latitude of the winter hemisphere neutral point. In both of 
these studies, however, the winter hemisphere neutral point moved outward 
and the summer hemisphere point moved inward, by approximately the same 


relative amount indicated in the present models. The substantially lower 
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latitude of the neutral points exhibited by the present models is probably 
another effect of the distributed currents within the magnetosphere, an 
effect which has been ignored in the various theoretically-derived models. 

Figure 4 shows field line projections in the dawn-dusk meridian for 
the quiet and disturbed models and tilt angles of 0^ and 30^. Substantial 
shifts in the high-latitude conjugate points occur for T = 30^- In general, 
the conjugate point in the winter hemisphere moves in longitude towards 
the midnight meridian. The line emerging at +68^ latitude, ±90^ solar 
magnetic longitude, for example, re-enters the earth at a longitude of 
±96.7^ for the SQ model and ±108^ for the SD model- There is very little 
shift in conjugate latitude, however , 

Figure 5 shows the field lines as viewed from above and projected 
into the equatorial plane for the quiet and disturbed models. As with 
previous field models, the field lines tend to be swept back into the 
tail, the effect increasing rapidly with increasing latitude. This effect 
is substantially greater at any given latitude than that predicted by the 
theoretical model of Mead (1964) ; compare his Figure 8) , but is consistent 
with the experimental results presented by Fairfield (1968) . It is likely 
that currents within the magnetosphere contribute substantially to this 
enhanced sweeping effect- 

Contours of the scalar difference Ab=|b| , 

^ 'model dipole' 

together with the field-line topology in the noon-midnight meridian plane, 
are shown in Figure 6 for the quiet and disturbed models and for T - 0^ 
and 30^. These may be compared with the experimentally-determined Ab 
contours given by Sugiura (l972a, Figure 1) and Sugiura et al . (1971, 
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Figures 1, 3), except that those studies did not attempt to determine 
the effect of varying tilt angle. Four general regions can be identified: 
(1) a +AB region of about 20-30Y near the subsolar point; (2) a -AB region 
of 20-30 y near the neutral points on the dayside, agreeing generally with 
the results of Fairfield and Mess (1972) ; (3) a high-latitude +AB region 

of 30-50Y on the nightside; and (4) a tongue-like -AB region of 10-25 y 
near the equator on the nightside. The general pattern agrees quite well 
with the results of Sugiura et al . except that they find a more intense 
and highly localized -AB region near the magnetic equator inside 5 
on both the dayside and the nightside, reaching -35 y even during quiet 
times. This featxare is not evident in our models because: (1) no data 

was used in our study inside 4 Rg, and most of the data between 4 and 
6 Rg was taken at high latitudes; and (2) our simple quadratic models 
are not capable of describing an extremely localized region of intense 
-AB. The high-latitude +AB regions on the nightside, however, is 
reproduced quite well by our models. For large tilt angles such as 
T = 30°, the +Ab regions are much more intense in the summer hemisphere, 
reaching 60-70Y/ than in the winter hemisphere, where Ab is only about 
20-40y. The B contours given by these models are much more realistic 
than those given by the Mead-Williams model (see Sugiura et al ., 1971, 
Figure 5) . 

The Ab contours in the dawn-dusk meridian are shown in Figure 7 . 

A 10-ISy depression is present near the equator, together with a 15-30Y 
region of positive Ab at higher latitudes and over the poles, which 
presumably connects with the high-latitude +Ab region on the nightside. 
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These contours compare favorably with those of Sugiura et al . {1971, 
Figures 2, 4) except for the intense -AB regions very near the equator. 
Figure 8 shows the contours of constant total field magnitude in 
the equatorial plane, as calculated with the quiet and disturbed model 
for T = 0°. As expected, the contours, which would be circles centered 
at the earth for a dipole field, are displaced toward the sun, indicating 
that the field is more intense on the dayside. The broad minimum on the 
nightside (|b| < IOy) is seen to extend around almost to the dawn and 
dusk meridians. These contours may be compared with those given by 
Fairfield (1968, Figure 2) , which were derived primarily from Explorer 21 
(Imp 2) data, except near the subsolar point, where Explorer 12 data was 
used. One discrepancy between these figures is that the earlier work 
gives field strengths near the svibsolar point roughly 20Y higher than the 
present work. A possible explanation of this discrepancy is that the 
solar wind flux may have been unusually high when Explorer 12 was near the 
subsolar point in 1961, thus producing magnetopause locations ( Cahill and 
Patel , 1967) closer to the earth than normal ( Fairfield , 1971), and 
leading to the higher field strengths in Fairfield's contours. If this 
is so, the contours in Figure 9 are more typical in this respect. A 
second discrepancy is the curvature of the contours near the dawn and dusk 
meridians in Figure 8, a feature which does not appear in Fairfield's 
contours. This curvature is probably an artificial effect, related to the 
tendency of the Z— component to increase with increasing distance beyond 
about 14 Rg. A third discrepancy is the presence of a steeper nightside 
gradient near 7 in Fairfield's contours. This may be a real effect. 
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which cannot be reproduced with our simple quadratic models - 

Figur e 9 shows the locus of the last closed field lines in the 
equatorial plane on the dayside for each of the four models for T = 0^, 
together with the average magnetopause position as observed by Fairfield 
(1971) . The calculated position for the two disturbed models agrees 
reasonably well with the observed magnetopause position. The last closed 
field line in the two quiet models, however, is 1-2 Rg outside the average 
observed magnetopause. The observations summarized by Fairfield (1971) 
were made over a wide range of tilt angles, whereas the model calculations 
were made with zero tilt. Calculations made with T ^ 0^ were somewhat 
more difficult to interpret, since in some regions near the expected 
magnetopause a field line traced in one direction will intersect the 
earth, whereas in the other direction it is open (e,g., see the -76 and 
-78^ field lines of the MF73SQ model in Figure 3) , However, even if one 
defines the magnetopause as the locus of the outermost point in the solar 
magneto spheric equatorial plane where only one end of the field line 
intersects the earth, this locus falls earthward of the last closed field 
line for T = 0^ except for a limited region near the subsolar point. 
Therefore, if one were to compute an average of the equatorial distance 
of the last closed field line over a range of tilt angles, this average 
would fall inside the corresponding T = 0° curves in Figure 9 and improve 
the agreement with the observations. 

It should be stressed at this point that no constraints other than 
those of symmetry have been placed on the models to achieve the various 
results described in this section. No observed or calculated magnetopause 
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positions have been used to derive the models • The only input are the 
vector magnetic field measurements inside the magnetopause expressed as 
functions of Y, Z, T, and Kp. The position of the last closed field 
line and neutral points on the day side ^ the minimum field region in the 
tail, the general alignment of field lines in the tail with the solar 
wind direction, etc., all result from making least-squares fits to the 
quadratic expansions. It is reassuring that the positions of the last 
closed field lines on the dayside agree as well as they do with the 
observed magnetopause positions. The agreement gives additional confidence 
in the general validity of the models. 

Figure 10 shows ATS*^5 measurements by Ski liman (1973) of the average 
diurnal variation in the three components of the field during quiet times 
at different seasons “of the year. The satellite was in a synchronous 
equatorial orbit at a longitude of 105^ west of Greenwich, placing it at 
a geomagnetic latitude of 9-3°. Also shown are the calculated values of 
the three components, using the IGRF internal field plus the MF73Q model. 

The calculated north or H component agrees quite well with the data for 
all three seasons, except for the average magnitude. Since the H com- 
ponent was along the satellite spin axis, the zero level was not precisely 
known, and a different choice for zero level would bring theory and 
experiment into better agreement. The predicted variation of the D or east 
component agrees remarkably well with the measurements at all seasons. The 
shape of the predicted vertical or V component differs somewhat from the 
measurements, but the average level and separation by seasons are in general 
agreement. Thus the model seems to predict the general characteristics of 
the field at 6.6 rather well. 
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The distributed magnetospheric currents implied by the models 

are easily calculated by taking the curl of Ek 5 uations 7-9» yielding 

(VxS) = (2c. - b, ) Y (16) 

X 4 1 

(VxB)y = (aj^ - c^) + (a^ - 2c^) X + (2a^ - Crj) T 2 (17) 

(VxB) = (b, - 2a^) T Y (18) 

(The curl of the dipole field is, of course, zero.) The values of each 

of these coefficients for the four models is given in Table 4. VxB is 

in units of y/lO^ when the solar magnetic coordinates are in units of 

tens of earth radii and T is in units of tens of degrees. By Maxwell's 

equations the ctirrent density is given by J = Vx§/y^ where 
—7 2 

= 4irxl0~^ newtons/amp in rationalized mks units. Thus to obtain 
the current density in amp/m^ the right-hand-side of Equations 16-18 
should be multiplied by 1.25x10 

Note that for T = 0° Equations 16-17 imply that the currents 
parallel to the equatorial plane are independent of Z and therefore not 
restricted to the plasma sheet in the tail. This unrealistic feature 
results from the fact that the field models contain only symmetric terms 
through quadratic in the spatial variables, and therefore the currents 
contain only linear terms. 

The current system in the X-Y plane for T = 0° implied by the 

MF73D model is shown in Figure 11. These currents are roughly concentric 

about X=2.5R,Y=0, and the current density increases linearly with 
2 

distance from this point. The large currents in the tail have a magnitude 
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of approximately lOy/R^ = 1.2xl0~^ amp/m^. This is comparable to the 
current densities obtained by Olson (1973) . If the tail has a diameter 
of approximately 30 Rg at this point, a rough calculation yields a value 

7 

of about 10 amperes for the total current flowing across the tail 
between X = -5 and X = -15 Rg. 

There is no indication of any ring currents near the 
earth such as that indicated in the work of Sugiura (1972a, 1972b, 1973). 
In principle, one should be able to determine the magnitude and extent 
of such a ring current by taking the curl of experimentally-derived 
models such as these. However, much more sophisticated models would 
have to be used in order to provide high-order terms in the Vxb expression 
appropriate to the geometry of such a ring current. The present 
expressions, containing only linear terms, are limited to determining 
only the gross current patterns in the outer magnetosphere. 
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MODEL LIMITATIONS 


A number of factors limit the accuracy and usefulness of these 

models. First of all, no data beyond 17 R was used as input to define 

E 

the models. Therefore, the models should not be used to calculate fields^ 

beyond 17 R„. Otherwise, unreliable results are likely to be obtained. 

E 

Secondly, there are large gaps in the data coverage. As indicated 
in an earlier section, 70% of the four -dimensional buckets contain no 
data. This lack of coverage is particularly evident at very high latitudes 
(especially high southern latitudes) , at near- equatorial latitudes, and at 
distances less than 5 There is also a lack of coverage at certain 

longitudes during some seasons of the year. The models may be 
inaccurate in these regions, just as models of the internal field based 
on surface measurements are likely to be inaccurate over large ocean and 
polar cap regions where little or no data are available. 

Thirdly, the present models cannot represent asymmetries in the 
magnetosphere. Only terms which satisfy the symmetry requirements of 
Equations 1-6 have been included. The models may therefore be suppressing 
real asymmetries, shch as the several-degree dawn -dusk asymmetries in 
declination recently reported by Sugiura (1973 ) - This particular 
asymmetry can be seen in our experimental data, but is suppressed by 
the models. 

Fourth, insufficient terms have been included in the models to 
represent the effects of relatively localized current systems. This 
is particularly true in the case of the near-earth equatorial ring 
current and sheet currents in the tail. Attempts to add higher-order 



29 


terms (for example, a 33-coefficient cubic model) gave rise to 
undesirable characteristics in regions where there was no data. 

Better data coverage is necessary before significant extensions of 
this type can be made. 

Finally, the models can only represent an average magneto- 
spheric configuration. Models such as these cannot represent the 
rapid time variations in magneto spheric configuration which take 
place during substorms. Such an extension to include time variations 
would obviously be highly desirable, but is undoubtedly far in the 
future. 

Despite these limitations, we believe that these models repre- 
sent the average field in the region from 4-17 R more accurately 
than any other model presently available. They have the additional 
advantage of being extremely simple to use; the FORTRAN decks 
representing the models are extremely short and require insignificant 
amounts of computer time to make a field calculation. Although other 
models exist which include either the effect of distributed currents 
or the effect of changes in the tilt angle, none contain both, and 
we feel that both effects must be included in any realistic model 
of the outer magneto spheric field. 



APPENDIX: COORDINATE SYSTEMS 


Equations 7“9 require as inputs solar magnetic coordinates and 
tilt angle- In many practical applications, however, the coordinates 
which are readily available are geographic latitude and longitude, day 
of year, and either Universal Time (UT) or local time. In addition, 
one frequently wishes to use an internal model such as the IGRF (lAGA, 
1969) together with these external models, rather than a simple dipole 
model. One must be careful in such cases to use a consistent coordinate 
system. Our preference is to make all such calculations in a geomag- 
netic dipole coordinate system, with the Z— axis parallel to the magnetic 
dipole. Appropriate formulas to transform spatial coordinates from 
geographic to geomagnetic latitude and longitude and vice versa are 
given by Mead (1970). The declination of the sun 6^ can be obtained 
from ephemeris tables, or an algorithm such as Mead*s given in Russell 
(1971) ; however, where high accuracy is not needed, the following 
algorithms are useful: 

sin = sin e sin (Al) 

Ag ~ 360 (IDAY - 81.6)/365.25 (A2) 

where e = 23.44^ is the obliquity of the ecliptic, is the sun’s 
celestial longitude along the ecliptic, and IDAY is the day of the year, 
with January 1 being day 1. The declination given by these algorithms, 
which assume that the sun moves uniformly along the ecliptic, is accurate 


to about half a degree. 





The tilt angle equal to the geomagnetic latitude of the sun, 
is defined by 

sin T= cos 11*4^ sin 6 + sin 11.4^ cos 6 cos (GHA^ - 69.8^) (A3) 

s s ^ 

where GHA^ is the Greenwich hour angle of the sun, given by 

GHA^ = GST - a =; 15 * UT - 180° (A4) 

S s ' ' 

where GST is Greenwich Sidereal Time and ot is the sun’s right ascension. 

s 

The position of the dipole (11.4° colatitude, 69.8° west longitude) is 
that given by the IGRF model at epoch 1965.0 ( Mead , 1970). Solar magnetic 
latitude is the same as geomagnetic latitude. Solar magnetic longitude 
is the same as magnetic local time (MLT) ± 12 hours. It is equal to the 
geomagnetic longitude, measured eastward from the dipole meridian, plus 
the magnetic hour angle of the sun, measured westward from the dipole 
meridian- The sun’s magnetic hour angle MHA^ is defined by 

sin MHA = cos 6^ sin (GHA - 69-8°) / cos T (A5) 

^ ^ s 

cos MHAg = (cos 11.4° sin T - sin 6^) / (sin 11.4° cos T) 

Equations A1-A6 plus those given by Mead (1970) can be used to 
calculate the solar magnetic coordinates and tilt angle needed as input 
to Equations 7—9. The internal field for the IGRF model can be obtained 
by using the coefficients given by Mead (1970) or Stern (1971) for input 
in magnetic dipole coordinates. (Fortran decks which include all these 
calculations are available from the authors.) 


(A6) 



Calculations of seasonal and diuirnal variations of geomagnetic 
phenomena must take into proper consideration the coherent variation 
of the tilt angle with local time during a 24'^hour day. This coherence 
will depend upon the geomagnetic longitude of the observer. At a dipole 
longitude of 0^ (69.8^ west geographic longitude) the tilt angle will be 
a maximum at noon and a minimum at midnight. The maximum and minimum 
are reversed at 180^ geomagnetic longitude, and a different phase 
relationship will exist at other longitudes. Thus even if a dipole is 
used for the internal field, the calculated diurnal variations of such 
phenomena as conjugate point positions will depend on the geomagnetic 
longitude of the observer as well as on the latitude and season of the 
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Table 1. Summary of the magnetometer data obtained 
from the four Imp spacecraft. 


Satellite 

Time 

Interva 1 

Total 

Time 

No. of 
Orbits 

No. of 
Individual 
Measurements 

No. of 
0 . 5Rg 
Averages 

Explorer 33 

7/66 - 7/68 

356 hrs . 

34 

0.25 X 10® 

744 

Explorer 34 

5/67 - 2/69 

1166 

144 

1.54 X 10® 

3510 

Explorer 41 

6/69 - 9/71 

1920 

195 

2.60 X 10® 

5668 

Explorer 43 

3/71 - 3/72 

898 

78 

40.37 X 10® 

2694 


Totals 4340 hrs. 451 44.76 x 10® 12616 



Table 2. Values of the coefficients to be substituted into Equations 


7-9 for each of the four models. The units are such that the contribution 


from 

each term 

is in gammas if the 

cartesian solar 

magnetic 

coordinates 

X, Y, 

and Z are 

in units 

of tens of 

earth radii and the tilt angle T 

is in 

units of 

tens of degrees. 




Compo^ 

nent 

Coeffi- 

cient 

Term 

Super- 

quiet 

MF73SQ 

Kp=0,0+ 

Quiet 

MF73Q 

Kp<2 

Disturbed 

MF73D 

Kp -2 

Super- 

Disturbed 

MF73SD 

Kp^3 

AB^ 

^1 

Z 

17-93 

21-79 

33.16 

39.48 



XZ 

-5.79 

-7.03 

-6.39 

-2.91 


^3 

T 

2-98 

3.02 

4.30 

5-17 


34 

TX 

-2.57 

-2.99 

-3. 25 

-3.86 



TX^ 

-0.30 

-0.62 

-0.44 

-1.04 



Ty2 

-1.47 

- 1.22 

-1.27 

-1.29 


^7 

T2|2 

1-05 

0.95 

0.45 

-1.14 

ABy 


YZ 

- 10.11 

-11.84 

-16.54 

-19.10 


^2 

TY 

-1.98 

-2.57 

-3.08 

-3.50 


b3 

TXY 

0.09 

- 0 . 28 

0.22 

0.23 



1 

-9.41 

-11.96 

-19.88 

-22.90 


^2 

X 

15.07 

17.87 

20.23 

22.70 


C3 


13.16 

15.88 

22.72 

26-50 


C 4 

Y^ 

8.36 

9.77 

13.23 

15.54 


^5 


7.95 

9.43 

11.46 

11.00 



TZ 

4.55 

5.57 

6.33 

7.36 



TXZ 

0.51 

1.53 

0.67 

1.85 





N-a 

Table 3. Root-mean- square values of the measured vector 
difference field (A b| „„_/ the vector residual 

L. lllo 

field |ab| and the figure-of -merit F (see 
text for definition ) , for each of the four models 



rms 

hSires 

F 

MF73SQ 

17.4 

8.2 

0.47 

MF73Q 

20.8 

10.1 

0.49 

MF73D 

29.9 

15.9 

0.53 


MF73SD 


34.8 


18.8 


0.54 



y/ 

Table 4. Coefficients for VxB in Equations 16-18. 
See text for discussion of units. 


Coefficients 

MF73SQ 

MF73Q 

MF73D 

MP73SD 

1 

o 

^1 

26.83 

31.38 

43.00 

50.18 

^1 

C2 

2.86 

3.92 

12.93 

16.78 


2C3 

-32.11 

-38.79 

-51.83 

-55.91 

2 3. ^ ■” 

^7 

1.59 

0.37 

0.22 

-4.13 

^3 - 


-3.03 

-2.16 

-2.76 

-2.81 




FIGURE CAPTIONS 

Figure 1- Distribution of the data set in four-dimensional 
coordinate space- The top half shows the distribution at 6-8 Rg 
and the bottom half at 13-15 Rg. The symbols indicate the range of 
the tilt angle associated with each measurement- 

Figure 2- Field lines and contours of constant total field magnitude 
in the noon-midnight meridian for T = Field lines are plotted 

every 2° in latitude / beginning with 66^- Contours are given every lOy 
from lO-lOOy- North- south symmetry is built into the models. Field 
lines and contours are terminated beyond 20 R^. 

Figure 3. Same as Figure 2, but for T = 30^. The symmetry is 
such that the northern and southern hemispheres would be exactly 
reversed for T = -30^. 

Figure 4. Field-line projections in the dawn-dusk meridian 
for T = 0^ and T = 30^- The field lines tend to be swept back into 
the tail. The view is along the solar magnetic equator. 

Figure 5. Field lines from rings of constant latitude, as viewed 
from above and projected into the equator. Lines are shown for every 
two hours of local time. Due to the basic symmetry of the models, 
the southern portion of each closed field line exactly overlaps the 
northern portion, and many of the lines which appear open are actually 
closed- Field lines from two models are plotted together; the lines 
from the MF73D model are always swept back further towards the 


tail. 



Figure 6- Contours of constant AB = |®j total ” 1^1 dipole' 
superimposed on the noon-midnight field topology- Contours are shown 
in lOy intervals. 

Figure 7. Contours of constant AB in the dawn-dusk meridian for 
T = and T = 30^- 

Figure 8- Contours of constant total field magnitude in the 
equatorial plane for T = O^. Contours are given every lOy from 
lOy-lOOy, The apparent increase in field magnitude beyond 13 Rg 
is probably not a real effect* 

Figure 9. Contours of the last closed field line in the equatorial 
plane given by the four models for T = plus the average position 
as observed by Fairfield (1971)* 

Figure 10. Diurnal and seasonal variation of the field at 
synchronous orbit. On the left are measured quiet- time averages 
obtained by Sugiura (1973) on ATS-5, located at a longitude of -105^. 

On the right are the values predicted by the MF73Q model at equinox 
and summer and winter solstice. The base line for the three components 
is that given by the IGRF internal model. 

Figure 11. Curl B in the equatorial plane for T = 0° (MF73D 
model). To get current densities in amp/m^z multiply the curl in 




























